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02CH1 1109,  a  graduate  assistantship  from  the  Department  of  Pharmaceutical  and  Biomedical 
Sciences  at  the  University  of  Georgia  and  a  subcontract  through  the  Henry  M.  Jackson 
Foundation  for  the  Advancement  of  Military  Medicine  (HJF)  from  the  Air  Force  Research 
Faboratory  (711  HPW/RHPB),  Wright-Patterson  Air  Force  Base  (WPAFB),  OH.  The  source  of 
funding  provided  by  RHPB  was  the  Environmental  Restoration  Account  (ERA)  project,  Military 
Site  Specific  Risk  Management,  Project  Number  CNBC20070001 .  Contract  oversight  for  HJF 
(contract  #  FA8650-05-2-6518)  was  provided  by  Mr.  Mark  Hoffman  of  71 1  HPW/RHP  at 
WPAFB,  OH.  The  authors  would  like  to  thank  Dr.  Melvin  E.  Andersen  for  his  technical 
assistance  with  model  development. 


SUMMARY 

Dichloroacetic  acid  (DCA)  is  a  common  disinfection  byproduct  in  surface  waters  and  a  probable 
minor  metabolite  of  triehloroethylene.  DCA  liver  eareinogenicity  has  been  demonstrated  in 
rodents  but  epidemiologieal  evidenee  in  humans  is  not  available.  High  doses  of  DCA  (10-50 
mg/kg)  are  used  to  treat  metabolie  aeidosis.  Biotransformation  of  DCA  by  glutathione 
transferase  zeta  (GSTzeta)  in  the  liver  is  the  major  elimination  pathway  in  humans.  GSTzeta  is 
inaetivated  by  DCA,  leading  to  slower  systemie  elearanee  and  nonlinear  pharmaeokineties  after 
multiple  doses. 

A  physiologieally-based  pharmaeokinetie  (PBPK)  model  was  developed  to  quantitatively 
describe  DCA  biotransformation  and  kinetics  in  humans  administered  DCA  by  intravenous 
infusion  and  oral  ingestion.  GSTzeta  metabolism  was  deseribed  using  a  Miehaelis-Menten 
equation  eoupled  with  rate  eonstants  to  aeeount  for  normal  GSTzeta  synthesis,  degradation  and 
irreversible  eovalent  binding  and  inhibition  by  the  glutathione-bound-DCA  intermediate.  The 
human  DCA  PBPK  model  adequately  predieted  the  DCA  plasma  kinetics  over  a  20,000  fold 
range  in  administered  doses.  Apparent  inhibition  of  GSTzeta  mediated  metabolism  of  DCA  was 
minimal  for  low  doses  of  DCA  (pg/kg/day),  but  was  signifieant  for  therapeutie  doses  of  DCA. 
Plasma  protein  binding  of  DCA  was  assumed  to  be  an  important  faetor  influeneing  the  kinetics 
of  low  doses  of  DCA  (pg/kg/day). 

Human  equivalent  doses  (HEDs)  were  ealeulated  for  a  1 0%  inerease  in  miee  hepatie  liver  caneer 
(2.1  mg/kg/day).  The  HEDs  for  the  dosimetries  (area  under  the  curve  (AUC)  for  total  and  free 
DCA  in  plasma,  AUC  of  DCA  in  liver  and  amount  of  DCA  metabolized)  per  day  were  0.02,  0.1, 
0.1  and  1.0  mg/kg/day,  respectively. 
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INTRODUCTION 

Dichloroacetic  acid  (DCA)  is  an  environmental  eontaminant  produeed  primarily  by  disinfeetion 
of  water  with  ehlorine.  The  reported  median  eoneentration  of  DCA  in  surfaee  water  systems  is 
15  pg/L  with  a  maximum  eoneentration  of  74  pg/L  (Boorman  et  al,  1999).  Coneentrations  as 
high  as  133  pg/L  may  also  oeeur  (Uden  and  Miller,  1983).  DCA  is  also  a  probable  minor 
metabolite  of  the  eommon  solvent  and  groundwater  eontaminant,  triehloro ethylene  (TCE). 

While  DCA  has  been  measured  as  a  by-produet  of  TCE  exposure  (Earson  et  al,  1992a),  these 
data  have  been  questioned  as  an  artifaet  of  sample  preparation  as  DCA  is  only  measurable  when 
high  levels  of  triehloroaoetie  aeid  (TCA)  are  present.  Merdink  et  al  (1998)  were  unable  to 
deteet  DCA  in  B6C3E1  miee  from  treatment  with  TCE  or  TCA  even  after  pre-treatment  with 
DCA  (2  g/E  in  drinking  water  for  2  weeks).  DCA  is  elassified  as  a  possible  human  eareinogen 
based  on  suffieient  evidenee  of  oaroinogenieity  in  rats  and  miee  (U.S.  EPA,  2003). 

Clinieally,  DCA  has  been  used  to  treat  metabolie  disorders  sueh  as  laetie  aeidosis  and  diabetes 
mellitus  due  to  its  ability  to  reduee  eireulating  glueose,  laetate  and  pyruvate.  Pharmaeologieal 
doses  of  DCA  range  from  10  to  50  mg/kg/day  (Barshop  et  al,  2004;  Jia  et  al,  2006).  Adverse 
effeets  from  therapeutie  use  inelude  mild  liver  dysfunetion,  transient  eentral  neuropathy, 
peripheral  neuropathy  and  hypoealeemia.  The  elinieal  effeets  are  generally  reversible  after 
withdrawal  of  treatments  (Staepoole  et  al,  1998b;  U.S.  EPA,  2003;  Mori  et  al,  2004). 

DCA  is  readily  absorbed  by  the  gastrointestinal  traet  as  only  1%  of  the  total  dose  is  exereted  in 
urine  after  a  single  oral  administration  of  50  mg/kg  in  humans.  The  plasma  half-life  of  DCA  in 
healthy  humans  is  approximately  1  hour  and  inereases  somewhat  with  repeated  dosing  (Curry  et 
al,  1991;  Staepoole  et  al,  1998a).  The  metabolism  of  DCA  involves  oxidative  deehlorination  to 
form  glyoxylate,  whieh  is  further  oxidized  to  oxalate,  earbon  dioxide  or  ineorporated  into  amino 
aeids  nueleophiles  and  other  eellular  moleeules  (Tong  et  al,  1998a;  U.S.  EPA,  2003;  Anderson 
et  al,  2004).  GSTzeta  is  the  only  enzyme  identified  in  the  biotransformation  of  DCA  (Tong  et 
al,  1998b;  Ammini  et  al,  2003;  Eim  et  al,  2004).  DCA  has  been  shown  to  inhibit  GSTzeta, 
through  eovalent  modifieation  of  the  enzyme,  in  both  in  vivo  animal  studies  and  in  vitro  studies 
(Anderson  et  al,  1999,  2002;  Sehultz  et  al,  2002).  GSTzeta  (GSTZl-1),  also  known  as 
maleylaeetoaeetate  isomerase  (MAAI),  is  an  essential  enzyme  in  phenylalanine/tyrosine 
eatabolism  pathway.  Disruption  of  tyrosine  eatabolism  by  eompetitive  inhibition  was  proposed 
as  a  possible  meehanism  for  DCA  toxieity  (Cornett  et  al,  1999;  Sehultz  et  al,  2002;  Ammini  et 
al,  2003;  Eantum  et  al,  2003).  Pour  funetional  allelie  variants  of  GSTZ  1-1  (GSTZ  la-la  and 
GSTZ  Id- Id)  have  been  identified  in  humans.  These  GSTZ  allelie  variants  have  different 
metabolie  eapaeities  and  inhibition  oharaeteristies  towards  DCA  in  vitro  (Tzeng  et  al,  2000; 
Blaekbum  et  al,  2001). 

The  first  rodent  physiologieally  based  pharmaeokinetie  (PBPK)  models  for  DCA  were  ereated  to 
traek  traee  amounts  of  DCA,  presumably  formed  as  a  metabolite  of  trichloroethylene  (Abbas  and 
Pisher,  1997;  Greenberg  et  al,  1999;  Clewell  et  al,  2000).  In  addition.  Barton  et  al  (1999) 
ereated  a  mouse  DCA  PBPK  model  to  understand  the  relationship  between  hepatie  eaneer 
ineidenee  and  DCA  hepatie  dosimetry  for  orally  ingested  DCA.  The  inhibitory  effeet  of  DCA  on 
its  own  metabolism  was  not  aeeounted  for  in  PBPK  models  until  Keys  et  al  (2004)  developed 
rodent  DCA  PBPK  models.  The  rodent  PBPK  models  of  Keys  et  al  (2004)  were  used  to 
evaluate  the  impaet  of  redueed  hepatie  metabolism  (by  suieide  inhibition  of  GSTzeta)  on  DCA 
blood  time  eourse  kineties. 

The  development  of  a  human  PBPK  model  for  DCA  was  possible,  in  part,  beeause  several 
pharmaeokinetie  studies  of  administered  DCA  have  been  reported  (Eukas  et  al,  1980;  Wells  et 
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al,  1980;  Chu,  1987;  Curry  et  al,  1991;  Fox  et  al,  1996;  Shangraw  and  Fisher,  1999;  Jia  et  al, 
2006;  Schultz  and  Shangraw,  2006).  In  the  present  paper  we  report  on  the  development  of  a 
human  PBPK  model  for  DC  A  formulated  to  account  for  reduced  hepatic  metabolism  of  DCA  via 
suicide  inhibition  of  the  GSTzeta  enzyme.  The  human  equivalent  doses  (HEDs)  for  a  10% 
increase  in  mouse  hepatic  liver  cancer  after  a  2  year  drinking  water  exposure  to  DCA  (2.1 
mg/kg/day)  were  then  estimated  using  the  human  DCA  PBPK  model. 


MATERIALS  AND  METHODS 

DCA  Human  Kinetic  Studies  Used  for  PBPK  Modeling 

Therapeutic  DCA  Dosing  Studies 

Curry  et  al.  (1991)  intravenously  infused  4  male  (average  weight  =  77  kg)  and  4  female  (average 
weight  =  54  kg)  healthy  subjects  over  a  30  minute  period  with  50  mg/kg  DCA.  Blood  samples 
were  collected  at  0.125  hour  (7.5  minutes)  intervals  for  the  first  hour  post-dosing,  at  0.25  hour 
(15  minutes)  intervals  for  the  second  hour,  at  0.5  hour  (30  minutes)  intervals  for  the  third  to  sixth 
hours  post-dosing,  and  then  at  1  hour  intervals  until  12  hours  post-dosing.  These  particular 
kinetic  studies  and  other  kinetic  studies  are  reported  in  a  dissertation  by  Chu  (1987).  Chu  (1987) 
collected  urine  from  each  of  these  8  subjects  at  various  intervals  up  to  12  hours  post-dosing  and 
calculated  the  cumulative  amount  of  DCA  excreted  and  DCA  urinary  clearance  rates.  In  other 
studies,  Chu  (1987)  reported  DCA  plasma  kinetics  after  a  single  oral  dose  of  50  mg/kg  DCA 
given  in  gelatin  capsules  to  4  male  subjects.  The  blood  sampling  schedule  was  similar  to  the 
intravenous  (IV)  experiments.  Chu  (1987)  also  administered  5  repeated  30  minute  IV  doses  of 
DCA  two  hours  apart  to  adult  volunteers.  Four  males  and  one  female  received  25  mg/kg  of 
DCA,  and  one  male  and  three  females  received  repeated  doses  of  50  mg/kg  DCA.  Blood 
samples  were  collected  at  30  minute  intervals  for  the  first  8.5-9  hours  after  the  first  infusion  and 
subsequently  every  3-4  hours  for  a  24  hour  period. 

Three  published  kinetic  studies  with  therapeutic  doses  of  DCA  were  used  to  test  the  ability  of  the 
model  to  predict  plasma  kinetics.  Lukas  et  al.  (1980)  intravenously  administered  either  10  or  20 
mg/kg  of  DCA  in  saline  to  two  adult  subjects  over  a  20  minute  period.  Blood  samples  were 
drawn  every  hour  for  the  first  4  hours,  then  every  2  hours  for  the  next  8  hours.  Shangraw  and 
Fisher  (1999)  treated  cirrhosis  patients  with  30  minute  IV  infusions  of  35  mg/kg  DCA.  Five 
healthy  men  and  one  healthy  woman  were  included  in  the  control  group.  In  these  individuals, 
blood  samples  were  drawn  at  0,  0.25,  0.5,  1,  2,  4,  6,  10,  14,  20  and  24  hours  post-dosing.  In  a 
somewhat  complicated  dosing  schedule.  Fox  et  al.  (1996)  administered  two  30  minute  IV 
infusions  of  DCA  at  8  hour  intervals  to  9  subjects  in  each  dose  group.  Doses  were  (1  St  ^  2"d 
infusion):  30+15,  60+30  or  100+50  mg/kg.  Blood  samples  were  collected  at  several 
predetermined  times  over  a  14  hour  period. 

Sub-therapeutic  DCA  Dosing  Studies 

Schultz  and  Shangraw  (2006)  administered  an  oral  dose  of  2  mg/kg  of  '^Ci-DCA  in  0.5  L  of 
water  followed  5  minutes  later  by  an  IV  dose  of  0.3  mg/kg  Ci-DCA  to  16  healthy  subjects. 
Blood  samples  were  collected  at  5  minute  intervals  for  the  first  0.5  hour,  10  minute  intervals  for 
the  remainder  of  the  first  hour,  and  at  1.25,  1.50  and  2.0  hours  post-dosing.  These  subjects  were 
then  kept  on  a  daily  oral  bolus  dose  of  0.02  mg/kg  of  ^^Ci-DCA  for  14  days  to  mimic  drinking 
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water  contaminated  with  DCA.  At  the  end  of  the  last  day  of  treatment  (Day  15),  the  subjects 
were  again  given  an  oral  bolus  dose  of  2  mg/kg  of  '^Ci-DCA  followed  by  an  fV  dose  of  0.3 
mg/kg  ‘^Ci-DCA. 

Jia  et  al.  (2006)  administered  2.5  |ag/kg  1,2-'^C-DCA  daily  to  20  healthy  adult  volunteers  (8 
male  and  12  female)  either  by  oral  or  IV  administration.  In  one  study,  15  of  these  volunteers  (6 
male  and  9  female)  received  orally  administered  for  either  5  or  15  consecutive  days  1,2-^^C- 
DCA  in  200  mL  distilled  water  after  an  overnight  fast.  Blood  was  sampled  on  Days  1,  5  and  15 
to  determine  1,2-^^C-DCA  plasma  kinetics.  In  the  2"^^  study,  16  volunteers  (7  male  and  9  female) 
received  daily  IV  doses  of  1,2-^^C-DCA  for  up  to  5  days.  Intravenous  1,2-^^C-DCA  (2.5  pg/kg) 
was  infused  at  a  constant  rate  in  saline  over  10  minutes  into  a  forearm  vein  after  an  overnight 
fast.  Blood  was  sampled  on  Days  1  and  5  to  determine  1,2-  C-DCA  plasma  kinetics.  On  each 
blood  sampling  day,  2  mL  of  venous  whole  blood  was  withdrawn  10  minutes  before 
administration  of  1,2-^^C-DCA,  and  then  0,  5,  10,  20,  30  minutes  and  1,  2,  3,  4,  6,  8,  12  and  24 
hours  after  administration  of  1,2-  C-DCA.  Urine  was  also  collected  at  12  and  24  hours  during 
each  kinetic  study.  Selected  volunteers  participated  in  both  the  oral  and  IV  dosing  study.  These 
1 1  volunteers  (5  male  and  6  female)  had  a  lapse  in  treatment  of  at  least  30  days  between  the  2 
routes  of  exposures. 

Human  PBPK  Model  for  DCA 

All  model  code  was  written  in  ASCL  (acslXtreme  V  2. 0.1. 2,  AEgis  Technologies  Group  Inc., 
Huntsville,  AL).  The  model  compartments  include  plasma,  liver,  kidney,  slowly  perfused  and 
rapidly  perfused  tissues  (Fig.  1).  Blood  flows  to  these  compartments  were  described  using 
venous  equilibration  equations.  Intravenous  dosing  was  described  as  an  infusion  rate  directly 
into  the  mixed  venous  blood  supply.  A  two-compartment  modeling  approach  (Abbas  and  Fisher, 
1997)  was  used  to  describe  oral  ingestion  of  DCA.  The  hepatic  metabolism  of  DCA  by  GSTzeta 
was  described  with  a  Michaelis-Menten  equation  modified  to  account  for  suicide  inhibition 
(Lilly  et  al,  1998).  Minor  elimination  of  DCA  through  the  urine  was  described  using  a  first- 
order  rate  constant.  Two  linked  PBPK  models  for  DCA  with  identical  structure  and  model 
parameter  values  were  used  to  describe  the  competitive  metabolism  of  two  isotopes  of  DCA  after 
co-administration  (Schultz  and  Shangraw,  2006).  Competitive  inhibition  of  metabolism  for 
chemical  mixtures  has  been  successfully  incorporated  into  PBPK  models  for  chemicals  that  are 
metabolized  by  the  same  enzymes  (Campbell  and  Fisher,  2007).  Suicide  inhibition  of  GSTzeta 
was  also  retained  in  both  models  where  co-administration  of  DCA  isotopes  was  simulated. 


4 


12C-DCA 


13C-DCA 


Figure  1.  Schematic  of  DCA  PBPK  model  for  humans. 

Each  isotope  was  descrihed  separately  to  allow  for  simulation  of  the  data  reported  hy 
Schultz  and  Shangraw  (2006),  Metabolic  interaction  was  added  to  the  liver  compartment 
to  maintain  the  separate  description  of  the  isotopes. 


The  metabolism  of  DCA  under  co-exposure  conditions  is  presented  in  Equations  1  and  2.  The 
amount  of  each  isotope  DCA  in  this  example)  that  is  metabolized  with  respect  to  time 
(dAm/dt,  mg/hour)  in  the  liver  is  given  by: 


RAM^^ 


V  X  Cvl" 

m  _ _  max  / _ 

dt  ~  K„x{l  +  Cv'^/KJ  +  Cv'^ 


(1) 


Assuming  no  isotope  effect,  the  Km  values  for  each  isotope  are  the  same,  thus  Equation  1  for 


DCA  can  be  rearranged  as: 


RAM^^  = 


K^+Cv]^  +  Cvf 


Similarly,  the  rate  of  metabolism  of  ^'^C  DCA  is  described  as: 

RAM'^  =  ^ 

dt  K^+Cvf+Cvf 


(2) 

(3) 


where  the  Michaelis-Menten  equation  represents  metabolic  conversion  by  the  GSTzeta  enzyme; 
Km  is  the  Michaelis-Menten  affinity  constant  (mg/E),  Vmax  (mg/hour)  is  the  current  metabolic 
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capacity  (inhibitable),  Cvf  and  Cvf  are  tbe  ^^C-DCA  and  ^^C-DCA  concentrations  in  tbe 
venous  blood  leaving  tbe  liver,  respectively. 

Tbe  proposed  meebanism  of  DCA-mediated  inaetivation  of  GSTzeta  involves  tbe  eovalent 
binding  of  tbe  glutatbione-bound-DCA  intermediate  metabolic  product  with  a  nucleopbilic  site 
on  GSTzeta.  This  covalent  modifieation  results  in  dose-dependent  inactivation  of  GSTzeta 
(Anderson  et  al,  2002).  In  tbe  present  human  model,  the  inactivation  reaction  rate  was 
mathematieally  deseribed  as  a  seeond  order  proeess  (bimoleeular  rate  eonstant,  kd,  1/mg),  whieh 
is  dependent  upon  the  intermediate  reaetive  metabolie  produet(s)  and  free  GSTzeta 
eoneentrations  in  the  liver  (Lilly  et  al,  1998).  The  rate  of  change  in  the  initial  value  of  maximal 
veloeity  of  metabolism  (Vmaxo)  with  respect  to  time  (dVmax/dt,  mg/h  )  is  shown  in  Equation  4. 

=  +  -k,  xV  ]  (4) 

T  .  d  V  -wy  '  X  jy  ^  ^  S  dc  ITiaX  J  X  Z 

at  +  Cvj  +  Cv, 

ml  ml 

2 

The  resynthesis  (kg,  mg/h  )  and  the  natural  degradation  (kde,  1/hour)  of  GSTzeta  (Equation  4) 
were  deseribed  in  a  similar  fashion  to  the  Keys  et  al.  (2004)  PBPK  model  for  rodents. 

Integration  of  Equation  4  provided  the  value  of  V^axas  a  function  of  inhibition  of  GSTzeta  using 
starting  maximal  veloeity  in  the  absenee  of  DCA  (Vmaxo). 
r‘dV 

Y  =  SS2L  + JZ  (-S') 

^  max  Jq  ^  maxO 

To  maintain  equilibrium  of  the  enzyme  at  steady  state,  k5=  kde  x  VmaxO-  The  pereentage  of 
reduction  of  the  GSTzeta  metabolie  eapaeity  (PVmax)  was  ealculated  as: 

Y 

PY^^=^^xm%  (6) 

IlldA  YP  ^ 

^  n 


where  Vmax(t)  is  the  maximum  metabolism  rate  at  time  t. 

This  PBPK  model  aeeounts  for  low-eapaeity,  high-affinity  binding  of  DCA  to  plasma  proteins. 
DCA  has  been  shown  to  be  weakly  bound  to  rat  plasma  protein  (e.g.,  6%,  Sehultz  et  al,  1999). 
Chu  (1987)  reported  that  23%  of  DCA  is  bound  in  human  plasma.  To  deseribe  the  rate  of 
binding  of  DCA  to  plasma  proteins  in  this  model,  the  methods  of  Clewell  et  al.  (2003)  were 
used,  where  the  binding  (assoeiation)  of  DCA  to  plasma  proteins  was  deseribed  using  Miehaelis- 
Menten  kineties,  and  the  rate  of  dissociation  (unbinding)  was  assumed  to  be  a  first-order  process. 
Physiological  parameters  used  in  the  human  model  are  summarized  in  Table  1.  Average  body 
weight  for  eaeh  pharmaeokinetie  study  was  used  when  reported.  Gender  differenees  in  DCA 
kineties  were  evaluated  in  seleeted  studies.  In  the  Shultz  and  Shangraw  (2006)  study,  DCA 
kinetic  data  were  reported  for  individual  subjects.  In  this  case,  physiologieal  model  parameters 
were  developed  for  males  and  females  separately  (Table  1)  and  also  for  both  sexes  eombined 
through  averaging  body  weight  and  blood  flows  to  the  liver  and  kidney,  whieh  are  slightly 
different  between  the  sexes.  No  appreeiable  gender  differenees  were  observed. 

Additional  differential  equations  used  in  the  human  model  and  a  list  of  parameter  name 
explanations  are  included  in  Appendix  A.  Appendix  B  eontains  the  model  eode  written  in 
aeslXtreme. 
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Table  1.  Physiological  model  parameter  values  used  in  human  PBPK  model  for  DCA 


Parameters 

Symbol 

Human 

Source 

Body  weight(kg) 

BW 

-70.0 

Subject  -specific  when  provided 

Cardiac  output 
(  l/h/kg^'^) 

QCC 

16.5 
15.87(7 
17.73  ? 

Gender-specific  when  provided. 

(Brown  etal.,  1997;  Fisher  etal.,  1998) 

Blood  flows  (  %  of  cardiac  output) 

Liver 

QLC 

26.5 

25.0(7 

27.0  ? 

Gender-specific  when  provided. 

(Brown  etal.,  1997;  Fisher  etal.,  1998) 

Kidney 

QKC 

17.5 

19.0  (7 

17.0  ? 

Gender-specific  when  provided. 

(Brown  et  al.,  1997;  Fisher  et  al.,  1998) 

Rapidly  perfused 

QRC 

32.0 

(Brown  et  al.,  1997;  Fisher  et  al.,  1998) 

Slowly  perfused 

QSC 

24.0 

(Brown  etal.,  1997;  Fisher  etal.,  1998) 

Tissue  Volumes 

(  %  of  body  weight) 

Plasma 

VP  LAC 

4.4 

(Merrill  etal.,  2005) 

Liver 

VLC 

2.6 

(Brown  et  al.,  1997) 

Kidney 

VKC 

0.44 

(Brown  et  al.,  1997) 

Rapidly  perfused 

VRC 

9.86 

(Brown  et  al.,  1997) 

Slowly  perfused 

VSC 

74.7 

(Brown  et  al.,  1997) 

Note:  h=hour 


Calibration  of  Human  PBPK  Model  for  DCA 

DCA  is  highly  hydrophilic  with  an  octanol/water  partition  coefficient  (logPo/w)  of  0.92. 
Chemieal-specific  parameters  for  DCA  are  listed  in  Table  2.  DCA  tissue/blood  partition 
eoefficient  (PC)  values  determined  for  mice  (Abbas  and  Fisher,  1997)  were  used  for  humans, 
except  where  noted.  A  small  urinary  clearanee  constant  (Cfc,  L/hour-kg)  value  was  ealeulated 
by  allometrieally  scaling  the  clearance  rates  reported  in  Curry  et  al.  (1991).  Cfc  was  fixed  prior 
to  further  model  parameterization. 

The  following  sequence  was  used  to  estimate  metabolie  parameters  describing  suicide  inhibition 
of  DCA  metabolism.  The  value  of  the  Miehaelis-Menten  affinity  eonstant  (Km)  was  set  to  6 
mg/L,  which  was  measured  using  human  liver  eytosol  (Tong  et  al,  1998b).  The  initial  value  of 
the  inhibition  rate  constant  (kd)  was  set  to  0.001  mg'^;  as  kd  and  Vmaxc  are  highly  eorrelated, 
setting  kd  to  a  very  small  value  allowed  initial  estimation  of  Vmaxc  using  low  dose  data.  The 
degradation  rate  (kde)  was  initially  set  to  the  value  used  for  rodents  (0.00875  hour'^;  Keys  et  al, 
2004).  It  was  necessary  to  lower  the  slowly  perfused  (PS)  partition  coefficient  value  obtained  in 
mice  (0.37,  Abbas  and  Fisher,  1997)  to  0.1 1  to  improve  the  agreement  between  observation  and 
model  prediction  following  0.3  mg/kg  ^^C-DCA  IV  bolus  on  Day  0  (Schultz  and  Shangraw, 
2006).  The  initial  estimate  for  the  maximal  rate  of  DCA  metabolism  (Vmaxco)  was  then  obtained 
by  visual  inspeetion  using  the  same  data  set  of  Schultz  and  Shangraw  (2006).  This  low  dose  of 
DCA  was  considered  to  have  a  minimal  effeet  on  GSTzeta  activity.  The  metabolie  inhibition 
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constant,  kd,  was  then  fit  to  the  high  dose  DCA  plasma  eoneentration-time  eourse  data  following 
5-repeated  infusions  of  50  mg/kg  and  25  mg/kg  DCA  at  2  hour  intervals  in  the  study  of  Chu 
(1987).  Fitting  was  aceomplished  by  the  maximum  likelihood  method  with  all  other  parameters 
fixed.  Next,  the  value  of  Vmaxco  was  again  visually  adjusted  to  fit  the  eoneentration-time  eourse 
data  from  a  50  mg/kg  single  IV  infusion  (Curry  et  al,  1991). 

Diminished  systemie  elearanee  of  DCA  after  multiple  or  single  high  doses  of  DCA  indieated  that 
all  individuals  (n  =12)  needed  more  than  4  weeks  to  reeover  (Chu,  1987;  Curry  et  al,  1991), 
suggesting  a  very  slow  rate  of  reeovery  for  GSTzeta  aetivity.  The  GSTzeta  resynthesis  rate,  kde, 
was  set  to  a  value  that  would  provide  a  reeovery  rate  for  GSTzeta  enzymatie  aetivity  equal  to  8 
weeks  after  a  single  50  mg/kg  dose  of  DCA  (2  out  of  4  subjeets  recovered  by  8  weeks).  The 
values  of  the  oral  absorption  rate  eonstants  following  drinking  water  administration  were 
estimated  by  fitting  the  C-DCA  eoneentration-time  eourse  data  after  2  mg/kg  C-DCA  oral 
ingestion  in  500  ml  of  drinking  water  reported  in  the  study  of  Sehultz  and  Shangraw  (2006).  The 
transfer  rate  from  the  first  GI  eompartment  to  the  liver  (Kai)  was  set  to  0.01  hour'^;  the  values  of 
the  transfer  rats  from  the  first  GI  to  second  GI  compartment  (Kai)  and  seeond  GI  eompartment  to 
the  liver  (Kas)  were  estimated  visually.  To  deseribe  low  eapaeity  binding  of  DCA  to  plasma 
proteins,  the  parameters  for  protein  binding  (Bmax,  Kmb  and  Kunb)  were  estimated  by  visually 
fitting  the  kinetic  data  from  IV  and  oral  administration  of  2.5  pg/kg  1,2-*^C-DCA  on  Day  1,  5 
and  15  (Jia  et  al,  2006).  Plasma  protein  binding  was  important  for  deseribing  the  plasma 
kinetics  of  DCA  after  administration  of  2.5  pg/kg  of  DCA. 
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Table  2.  Chemical  specific  parameter  values  for  human  PBPK  model  for  DC  A 


Parameters 

Symbol 

Value 

Source 

Partition  coefficient 

LiverBlood 

PL 

1.08 

Abbas  and  Fisher,  1997 

Kidney:Blood 

PK 

0.74 

Abbas  and  Fisher,  1997 

Rapidly  perfused:Blood 

PR 

1.08 

Abbas  and  Fisher,  1997 

Slowly  perfused:Blood 

PS 

0.11 

Estimated  (Schultz  and  Shangraw, 
2006) 

Metabolic  and  clearance  parameters 

Maximal  rate  of  GSTzeta  velocity 

V  maxcO 

50 

Estimated  (Curry  et  a!.,  1991; 

(mg/h/kg°'"®) 

Schultz  and  Shangraw,  2006) 

Michaelic-Menten  affinity  constant 
(mg/I) 

Km 

6.0 

Tong  et  al.,  1998a 

Inhibition  constant  (  mg 

kd 

0.004 

Estimated  (Chu,  1987) 

Degradation  rate  constant  (  h'^) 

kde 

0.004 

Estimated  (Curry  ef  a/.,  1991) 

Urinary  excretion  (l/h-kg) 

Cfc 

0.7x10'^ 

Curry  et  al.,  1991 

Oral  absorption  parameters  following  drinking  water  ingestion 


Transfer  rate  from  the  C*  Gl 
compartment  to  the  liver  (  h'^) 

Ka1 

0.01 

Estimated  (Schultz  and  Shangraw, 
2006) 

T ransfer  rate  from  the  1  to  the 

2"^^  Gl  compartment  (  h'^) 

Ka2 

7.0 

Estimated  (Schultz  and  Shangraw, 
2006) 

Transfer  rate  from  the  2"'^  Gl 
compartment  to  the  liver  (  h'^) 

Ka3 

7.0 

Estimated  (Schultz  and  Shangraw, 
2006) 

Plasma  protein  binding  parameters 

Maximum  capacity  (mg) 

^max 

0.06 

Estimated  (Jia  et  al.,  2006) 

Affinity  constant  (mg/I) 

^mb 

0.001 

Estimated  (Jia  et  al.,  2006) 

Dissociation  constant  (h'^) 

^unb 

0.16 

Estimated  (Jia  et  al.,  2006) 

Note:  h=hour 


Individual  metabolic  capacities  were  estimated  for  the  16  subjects  reported  in  Schultz  and 
Shangraw  (2006)  to  assess  the  inter-individual  variation  as  inferred  by  differences  in  clearance 
kinetics  of  DCA  across  the  study  participants.  Gender  specific  values  of  physiological 
parameters  were  used  (Table  1)  but  not  considered  essential  because  gender  appears  to  have  little 
influence  on  DCA  kinetics  (Stacpoole  et  al,  1998a;  Schultz  and  Shangraw,  2006).  Vmaxco,  kd  and 
the  oral  absorption  parameters  (Ka2  and  Ka3)  for  each  subject  were  estimated  by  maximum 
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likelihood  methods.  First,  the  initial  estimations  of  VmaxcO,  kd,  Ka2  and  Ka3  were  fit 

12  1 3 

simultaneously  to  the  individual  C-DCA  and  C-DCA  eoneentration-time  eourse  data  of  eaeh 
subjeet.  Seeondly,  Ka2  and  Kas  were  fixed  at  the  initial  estimated  values  and  Vmaxco  and  kd  were 

re-optimized  by  htting  to  C-DCA  and  C-DCA  eoneentration-time  eourse  data.  Finally,  Ka2 

12 

and  Ka3  were  re-optimized  by  fitting  to  C-DCA  eoneentration-time  eourse  data  with  ah  other 
parameters  fixed. 


Predictive  Performance  Analysis 

The  ability  of  the  PBPK  model  to  prediet  the  pharmaeokineties  of  DCA  in  humans  for  eaeh  data 
set  was  evaluated  using  deseriptive  measures  that  eompare  prediction  with  observation 
(Gustafson  et  al,  2002).  The  predictive  performance  measures  were:  1)  ratio  of  the  PBPK  model 
predicted  area  under  the  concentration  curve  (AUC)  for  the  total  plasma  concentration  of  DCA 
divided  by  the  experimental  AUC  for  the  total  plasma  DCA  measurements  calculated  by  non- 
compartmental  analysis  (WinNonlin  V4.1,  Pharsight,  Mountain  View,  CA)  (AUC(o^t)  P/M 
Ratio,  Equation  7);  2)  the  performance  errors  (PEs,  Equation  8)  and  the  median  value  for  ah  the 
calculated  PEs  (MPE%,  Equation  9);  3)  the  median  value  for  the  absolute  value  of  ah  calculated 
PEs  (MAPE%,  Equation  10);  and  4)  the  root  mean  squared  PEs  (RMSPE%).  The  calculations 
are: 


AUC 


(O^t) 


AUC 

P/M  Ratio  =  predicted 


AUC. 


(V) 


( 0  — >  ? )  exp  erim  en  tal 


Performance  Error  (PE)  = 


C  -C 

measured  predicted  i/x/xn/ 

— - xl00% 


c 


predicted 


MPE%  =  Median  (PEi,  PEi,...  PEn)  (9) 
MAPE%  =  Median  (  |PEi|,  IPE2I,  ...|PEn|  ) 


RMSPE%  = 


(ll) 


(lO) 


(8) 


where  AUC(o^t)predicted  and  AUC(o^t)experimentai  arc  the  PBPK  model  predicted  and  experimental 
plasma  AUCs  until  the  last  measured  data  point,  respectively.  Cmeasured  and  Cpredicted  are  the 
measured  and  model  predicted  DCA  plasma  concentrations,  |PEi|  is  the  absolute  value  of  PEi  and 
i  =  I,  2....n,  representing  the  sample  size  in  each  data  set.  MPE%  is  a  measure  of  bias  and 
MAPE%  and  RMSPE%  are  measures  of  accuracy  (Gustafson  et  al,  2002). 


Sensitivity  Analysis  of  Parameters 

A  sensitivity  analysis  using  the  automated  process  provided  in  acslXtreme  was  carried  out  using 
AUC  for  DCA  in  plasma.  Briefly,  the  equation  used  to  calculate  the  sensitivity  coefficient  for 
each  parameter  is: 

Sensitivity  Coefficient  = ^  (12) 

(C-D)/D 

where  A  is  the  DCA  plasma  AUC  predicted  with  a  1%  increased  parameter  value,  B  is  the  DCA 
plasma  AUC  predicted  at  the  starting  parameter  value,  C  is  the  parameter  value  after  a  1% 
increase  and  D  is  the  original  parameter  value.  A  large  absolute  value  of  the  coefficient  indicates 
the  parameter  has  an  important  effect  on  the  model.  A  positive  value  indicates  that  the  model 
parameter  and  the  corresponding  output  are  positively  related  and  a  negative  value  indicates  they 
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are  inversely  related.  The  ealculation  was  performed  using  the  forward  differenee  algorithm 
with  delta  equal  to  0.01  and  normalized  to  both  response  variables  and  parameters.  Initial 
calculations  were  checked  by  hand  to  verify  the  algorithm  output. 

Cancer  Risk  Estimates 

The  theoretical  cancer  risk  for  ingestion  of  DC  A  was  calculated  using  a  previously  reported 
mouse  PBPK  model  for  DCA  (Keys  et  al,  2004)  and  the  current  human  PBPK  model  for  DCA 
assuming  a  70  kg  individual  ingests  2  L  of  water  each  day.  One  liter  is  ingested  over  a  5  minute 
period  every  12  hours.  A  point  of  departure  (POD)  for  a  10%  increase  in  cancer  risk  for  DCA  is 
reported  to  be  2.1  mg/kg/day  (estimated  as  the  95%  lower  confidence  limit  on  the  benchmark 
dose)  (U.S.  EPA.,  2003)  based  on  liver  carcinomas  in  B6C3Fi  mice  exposed  to  DCA  in  drinking 
water  (DeAngelo  et  al,  1999).  The  previously  published  DCA  mouse  model  (Keys  et  al,  2004) 
was  used  to  estimate  the  dosimetries:  1)  daily  AUC  of  free  plasma  DCA  (mg/L/day);  2)  daily 
AUC  of  total  plasma  DCA  (mg/L/day);  3)  daily  liver  AUC  (mg/L/day);  and  4)  the  daily  amount 
of  DCA  metabolized  per  kg  of  liver  (Ami,  mg/kg/day)  at  steady  state  following  consumption  of 
2.1  mg/kg/day  DCA  in  drinking  water.  Assuming  a  70  kg  human  ingests  1  L  every  12  hours,  the 
present  human  model  was  used  to  estimate  human  equivalent  doses  (HEDs)  corresponding  to  the 
4  dosimetries  in  mice. 

The  rationale  for  using  AUC  for  DCA  in  the  liver  and  the  amount  of  DCA  metabolized  in  the 
liver  is  related  to  a  lack  of  information  pertaining  to  the  specific  mechanism  of  action  for  DCA 
liver  carcinogenicity.  DCA,  reactive  by-products  from  metabolism,  or  a  combination  of  the 
parent  and  reactive  metabolites  may  be  responsible  for  the  resulting  liver  cancer.  These  two 
dosimetries  are  good  measures  of  internal  dose  because  liver  is  the  target  organ  and  measured 
liver  DCA  concentrations  were  available  with  which  to  validate.  However,  when  extrapolating 
to  humans,  measurements  of  DCA  in  liver  are  not  available,  but  measurements  of  DCA  are 
available  in  plasma.  The  plasma  concentration  represents  the  dose  available  for  delivery  to  the 
liver  and,  hence,  is  a  good  surrogate  dosimetric  for  DCA  liver  concentration. 

In  this  case  we  speculate  that  DCA  is  bound  to  plasma  proteins  at  low  capacity  and  high  affinity 
and  the  bound  DCA  is  potentially  not  available  for  uptake  into  the  liver.  As  there  is  no  direct 
supporting  experimental  evidence,  we  calculated  dose  metrics  for  total  and  free  DCA  plasma 
concentrations. 

RESULTS 

The  range  of  administered  IV  and  oral  ingestion  doses  of  DCA  described  with  this  PBPK  model 
was  a  remarkable  20,000  fold  (2.5  pg/kg  to  50  mg/kg)  with  studies  conducted  at  several 
laboratories.  Also  noteworthy,  one  set  of  model  parameter  values  (Tables  1  and  2)  was  used  for 
the  DCA  PBPK  model  to  describe  this  entire  data  set  with  the  exception  of  individual  subject 
optimized  kinetic  data  reported  by  Schultz  and  Shangraw  (2006). 

Perhaps  the  most  challenging  pharmacokinetic  data  to  describe  with  the  present  PBPK  model 
were  the  high  dose  DCA  studies  conducted  in  the  laboratory  of  Chu  (1987).  The  model 
predicted  systemic  clearance  of  DCA  over  a  6  hour  period  was  over-predicted  in  4  male  and  4 
female  subjects  (Fig.  2)  after  a  single  IV  infusion  of  50  mg/kg  (Curry  et  al,  1991)  and  modestly 
under-predicted  following  5 -repeated  administrations  of  DCA  (50  and  25  mg/kg)  at  2  hour 
intervals  over  a  10  hour  period  (Fig.  3).  The  cumulative  urinary  excretion  of  DCA  was 
somewhat  variable  (Fig.  4)  in  these  subjects,  accounting  for  less  than  1%  of  the  administrated 
dose.  The  PBPK  model  predicted  cumulative  amount  of  DCA  excreted  in  urine  is  in  good 
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agreement  with  observed  exeretion,  but  is  slightly  below  the  mean  eumulative  amount  of  urine 
observed  in  these  subjects  (Fig.  4). 


0  2  4  6 

Time  (hrs) 

Figure  2,  Model  predicted  ( — )  vs.  observed  (•)  mean  DCA  plasma  concentrations  for  4 
female  subjects,  and  model  predicted  ( — )  vs.  observed  (o)  mean  DCA  plasma 
concentrations  for  4  male  subjects  following  a  0.5  hour  IV  infusion  of  50  mg/kg  DCA.  The 
measured  values  were  originally  reported  by  Curry  et  al.  (1991). 
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DCA  in  Plasma  (mg/L) 


Figure  3,  Model  predicted  and  observed  mean  DCA  plasma  concentrations  for  24  hours 
following  5-repeated  IV  0,5  hour  infusions  at  2  hour  intervals.  The  concentration-time 
course  data  for  2  doses  are  presented:  25  mg/kg  (o,  n=5)  and  50  mg/kg  (•,  n=4).  The 
measured  data  were  originally  reported  by  Chu  (1987), 
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60 


50  A 


Time  (hrs) 

Figure  4,  Model  simulation  ( — )  of  the  cumulative  amount  of  DCA  excreted  in  the  urine  vs. 
the  average  cumulative  excretion  measured  in  8  subjects  following  a  0.5  hour  IV  infusion  of 
50  mg/kg  DCA.  Measured  values  were  originally  reported  hy  Chu  (1987).  Bars  represent 

standard  deviation. 


For  other  high  dose  studies  (Lukas  et  al,  1980),  with  2  individuals  eaeh  intravenously 
administered  10  or  20  mg/kg,  good  agreement  was  obtained  between  model  predietions  and 
observations  (Fig.  5).  In  another  group  of  6  individuals  administered  35  mg/kg  of  DCA 
(Shangraw  and  Fisher,  1999),  the  systemie  elearanee  of  DCA  was  slightly  over-  predicted  (Fig. 
5).  For  the  final  high  dose  simulation  study  (Fig.  6)  in  which  groups  of  individuals  were 
intravenously  administered  a  second  dose  of  DCA  8  hours  after  the  first  dose  (Fox  et  al,  1996), 
favorable  agreement  was  obtained  between  model  predictions  and  observations. 
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Figure  5.  Model  predicted  ( — )  and  observed  DC  A  plasma  concentrations  of  2  subjects 
following  10  mg/kg  (•jA)  and  2  subjects  following  20  mg/kg  (o,  A)  IV  infusions  of  DCA 
administered  over  20  minutes  (originally  reported  by  Lukas  et  aL,  1980),  Also  shown  are 
the  mean  plasma  concentrations  (x,  n=6)  following  a  0,5  hour  IV  infusion  of  35  mg/kg  DCA 

reported  by  Shangraw  and  Fisher  (1999), 
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Figure  6,  Model  predicted  ( — )  and  observed  mean  DCA  plasma  concentrations  after  2  0,5 
hour  IV  infusions  of  DCA  at  0  and  8  hours.  The  first  and  second  administered  doses  of 
DCA  were  30  +  15  mg/kg  (•,  n=9),  60  +  30  mg/kg  (o,  n=9)  and  100  +  50  mg/kg  (A,  n=9), 
respectively.  Experimental  data  were  originally  reported  by  Fox  et  al  (1996), 


Two  independent  elinieal  studies  were  used  for  low-dose  DCA  simulation  studies.  In  the  Sehultz 

12 

and  Shangraw  (2006)  study,  individuals  were  initially  dosed  orally  with  2  mg/kg  of  C-DCA 

1  -5 

followed  within  minutes  by  an  IV  dose  of  0.3  mg/kg  of  C-DCA.  After  14  days  of  daily 
treatment  with  a  low  dose  (0.02  mg/kg)  of  ^^C-DCA  in  drinking  water,  this  oral/IV  dosing 
sehedule  was  repeated  and  serial  blood  samples  were  eollected.  Fig.  7a  depiets  the  entire 
simulation  of  both  oral  and  IV  dosing  of  men  at  the  beginning  and  end  of  the  study  as  well  as  the 
daily  low-dose  treatment.  DCA  plasma  time  eourse  data  for  both  the  IV  and  oral  administration 
are  grouped  tightly.  In  panels  b  and  e  of  Fig.  7,  observed  and  model  predieted  plasma 
eoneentrations  are  shown  for  a  2  hour  period  after  oral  and  then  IV  administration  of  the  two 
isotopes  of  DCA  at  the  beginning  and  end  of  the  treatment  period.  Generally  speaking,  good 
agreement  was  observed  between  model  predieted  DCA  eoneentrations  and  measured  values. 

The  model  predicted  clearance  of  IV  administered  DCA  in  men  was  somewhat  slower  than 
observed  for  both  days  on  which  dosing  occurred. 

Similar  trends  between  observed  data  and  model  prediction  were  found  for  women  (Fig.  8). 

When  individual  DCA  clearance  kinetics  were  evaluated  by  adjusting  the  values  for  metabolism 
(Vmaxco,  kd)  and  oral  absorption  (Ka2  and  Kas)  (Table  3),  the  fitted  values  for  VmaxcO  had  the 
largest  range  of  estimated  values  (36.7  to  280  mg/hour/kg°'^^),  while  values  for  kd  ranged  from 
zero  (no  metabolic  inhibition)  to  0.016  mg'\  The  individual  values  of  Ka2  and  Ka3  are  relatively 
similar  and  the  average  values  are  about  8.0  hour’'.  The  simultaneously  optimized  values  for 
VmaxcO,  kd,  Ka2  and  Ka3  wcrc  very  close  to  the  final  values  obtained  after  the  re-optimization  was 
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conducted  separately  for  the  metabolic  parameters  (VmaxcO,  kd)  and  the  oral  absorption  parameters 
(Ka2  and  Kas).  This  was  expected  because  the  metabolie  parameters  influence  the  AUC  of  both 
C-DCA  and  C-DCA  while  the  oral  uptake  parameters  only  influence  the  shape  of  C-DCA 
eoneentration-time  profiles  (i.e.,  peak  eoneentration  (Cmax)  and  time  to  reaeh  the  peak 
coneentration,  Tmax). 

In  the  Jia  et  al.  (2006)  clinieal  study,  individuals  were  administered  a  very  low  dose  of  2.5 
pg/kg/day  for  5  days  by  IV  infusion  or  orally  for  15  days.  Blood  was  collected  on  speeified  days 
following  treatment.  These  low-dose  kinetic  data  demonstrated  the  requirement  for  plasma 
protein  binding  of  DCA  in  the  model  (Fig.  9a),  as  the  kinetic  behavior  eould  not  be  described 
without  assuming  plasma  protein  binding  (Fig.  9b).  In  the  other  higher-dose  DCA  studies,  the 
influence  of  this  low-capacity,  high-affinity  plasma  protein  binding  on  DCA  kineties  was  not 
readily  apparent.  The  binding  of  DCA  to  these  unspeeified  plasma  proteins  was  assumed  to  be 
reversible.  When  2.5  pg/kg/day  of  DCA  was  administered  orally,  peak  plasma  DCA 
concentrations  were  moderately  under-predicted  on  Days  1,  5  and  15  of  treatment  (Fig.  10). 
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Figure  7.  Model  predictions  of  ^^C-DCA  ( — )  and  ^^C-DCA  ( — )  vs.  observed  mean  ^^C- 
DCA  (A)  and  ^^C-DCA  plasma  concentrations  (•)  in  8  male  subjects.  Panel  (a)  includes  the 
entire  period  of  the  study  (15  days)  where  subjects  were  administered  an  initial  oral  dose  of 
2  mg/kg  ^^C-DCA  in  500  mL  of  water  followed  immediately  hy  an  IV  holus  of  0,3  mg/kg  of 
^^C-DCA,  The  subjects  were  maintained  on  a  daily  dose  of  0,02  mg/kg  ^^C-DCA  in  500  mL 
of  water  for  14  days  after  which  time  the  oral/IV  dosing  regime  was  repeated  on  Day  15, 
Panel  (b)  shows  the  initial  2  hours  following  the  first  2  mg/kg  oral  and  0,3  mg/kg  IV  doses 
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on  Day  0  and  panel  (c)  gives  the  initial  2  hours  following  the  second  2  mg/kg  oral  and  0,3 
mg/kg  IV  administration  on  Day  15.  The  data  were  originally  reported  in  Schultz  and 

Shangraw  (2006). 


Time  (hrs) 

Figure  8,  Model  predictions  of  ^^C-DCA  ( — )  and  ^^C-DCA  ( — )  vs.  observed  mean  ^^C- 
DCA  (A)  and  ^^C-DCA  plasma  concentrations  (•)  in  8  female  subjects.  Panel  (a)  gives  the 
initial  2  hours  following  the  first  2  mg/kg  oral  and  0,3  mg/kg  IV  administration  on  Day  0, 
and  panel  (b)  gives  the  initial  2  hours  following  the  second  2  mg/kg  oral  and  0,3  mg/kg  IV 
administration  on  Day  15  (see  Fig.  7  for  detailed  dosing  explanation).  Data  originally 

reported  in  Schultz  and  Shangraw  (2006), 
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Table  3,  Optimized  parameter  values  for  metabolic  constants  (Vmaxco  and  kd)  and  oral 
uptake  constants  (Ka2  and  Ka3)  using  individual  DCA  plasma  pharmacokinetic  data  from 

Shultz  and  Shangraw  (2006) 


Subject 

Gender 

Weight 

(kg) 

VfnaxcO 

(mq/h/kg'*-^®) 

kd 

(mg^) 

Ka2 

(h-^) 

Ka3 

(h-^) 

1 

(? 

90 

113.2 

0.001 

5.1 

5.1 

2 

? 

66 

71.1 

0.000 

5.3 

5.3 

3 

? 

74.5 

118.8 

0.011 

9.3 

9.1 

5 

(? 

66.6 

113.2 

0.002 

5.6 

2.9 

6 

9 

78.3 

252.2 

0.016 

6.9 

6.9 

7 

(9 

64.2 

89.3 

0.000 

4.6 

4.0 

8 

9 

59 

61.4 

0.005 

7.9 

7.9 

9 

(? 

79 

97.9 

0.005 

9.0 

9.0 

10 

9 

45 

225.8 

0.011 

7.7 

12.5 

11 

9 

69.5 

36.7 

0.006 

9.5 

9.5 

12 

(? 

78.6 

280.0 

0.010 

8.8 

8.9 

13 

(? 

73.5 

58.4 

0.000 

11.8 

12.0 

15 

9 

67.4 

60.0 

0.006 

17.7 

5.7 

16 

(9 

55.9 

83.1 

0.000 

8.0 

8.1 

17 

9 

67.8 

51.1 

0.000 

7.7 

13.3 

18 

(? 

73.2 

45.8 

0.000 

8.9 

8.9 

Mean 

69.3 

109.9 

0.005 

8.4 

8.1 

SD 

10.6 

75.7 

0.005 

3.1 

3.0 
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Figure  9,  Model  predicted  ( — )  and  observed  (•)  ^^C-DCA  concentrations  following  a  10 
minute  IV  infusion  of  2,5  pg/kg/day  for  5  days.  Panel  (a)  shows  the  final  model  with  the 
added  low-capacity  plasma  protein  binding  and  panel  (b)  shows  the  initial  model  without 
binding.  The  measured  data  originally  were  reported  by  Jia  et  al  (2006), 
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Time  (hrs) 

Figure  10,  Model  predicted  ( — )  and  observed  (A)  ^^C-DCA  concentrations  following 
ingestion  of  2,5  pg/kg/day  DCA  in  200  mL  of  water  for  15  days.  Measured  data  originally 

were  reported  by  Jia  et  al.  (2006), 


Prediction  Performance  Analysis 

In  addition  to  visual  inspection,  the  goodness-of-fit  of  a  model  prediction  to  each  data  set  was 
evaluated  with  the  AUC  P/M  Ratio,  MPE%,  MAPE%  and  RMSPE%  (Table  4).  The  AUC  P/M 
Ratio  values  ranged  from  0.4  to  2.9  for  23  pharmacokinetic  data  sets  with  a  mean  AEIC  P/M 
Ratio  value  of  1.08.  A  value  near  1.0  indicates  a  good  agreement  between  observation  and 
prediction.  The  values  of  MPE%  ranged  from  -50.6  to  202.2%.  Thirteen  of  the  23 
pharmacokinetic  data  sets  had  MPE%  values  that  were  negative  and  the  remaining  10 
pharmacokinetic  data  sets  were  positive.  The  even  distribution  of  MPE%  values  around  zero 
suggests  there  was  little  bias  in  the  PBPK  model  to  either  under-  or  over-  predict  the 
experimental  data  systematically.  MAPE%  and  RMSPE%  values  (Table  4)  for  the  23 
pharmacokinetic  data  sets  suggest  that  a  few  of  the  PBPK  model  predictions  (Jia  et  al,  2006; 
Shangraw  and  Eisher,  1999)  were  less  accurate,  as  indicated  by  large  MAPE%  and  RMSPE% 
values  relative  to  the  values  for  the  other  pharmacokinetic  data  sets.  These  performance  metrics 
appear  to  give  more  weight  to  the  data  points  with  smaller  values  because  the  calculations  are 
based  on  relative  percent  error  rather  than  absolute  percent  error. 

Sensitivity  Assay 

Since  the  sensitivity  analysis  is  both  dose  and  time  dependent,  this  analysis  was  applied 
using  dosimetries  that  reflect  integrated  measures  over  time,  such  as  AEIC  (Table  5). 

Normalized  sensitivity  coefficients  with  absolute  values  that  were  determined  to  be  most 
sensitive  (>0.5)  are  in  bold.  Cardiac  output  (QCC)  and  blood  flow  to  the  liver  (QEC)  were  found 
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to  be  sensitive  for  low  IV  doses  of  DCA  (0.3  mg/kg  and  2.5  |ag/kg),  while  the  metabolie 
parameters  (Vmaxco,  and  kd)  were  the  most  sensitive  parameters  at  therapeutic  doses  (10  to  50 
mg/kg/day  (Barshop  et  al,  2004;  Jia  et  al,  2006)).  The  protein  binding  parameters  (Bmax  and 
Kunb)  were  sensitive  model  parameters  only  for  the  lowest  DCA  dose  of  2.5  pg/kg/day. 
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Table  4.  Predictive  Performance  of  Model  Predictions  by  AUC(o_.t)  P/M  Ratio,  MPE%,  MAPE% 
and  RMSPE% 
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100+50  0.9  4.5  5.8  38.9 


Table  5,  Normalized  sensitivity  coefficients  of  selected  parameters  in  the  PBPK  model  for 
DCA,  based  on  model  predictions  of  DCA  plasma  AUC  in  various  exposure  scenarios 


Model 

Parameters 

Normalized  Sensitivity  Coefficients 

50mg/kg 

iv 

(single) 

AUCo^24h 

50mg/kg 

iv 

(2  h  interval) 

AUCo^24h 

0.3mg/kg 

iv 

(Day  0) 

AUCo^24h 

2mg/kg 
oral 
(Day  0) 
AUCo^24h 

2.5  |jg/kg-d 
iv 

AUCo^i20h 

2.5  |jg/kg-d 
oral 

AUCo^360h 

QCC 

-0.05 

-0.02 

-0.54 

-0.02 

-0.64 

-0.15 

QLC 

-0.09 

-0.03 

-0.50 

_ 

-0.57 

-0.12 

V  maxcO 

-0.88 

-0.53 

-0.50 

-1.05 

-0.33 

-0.72 

Km 

0.47 

0.41 

0.43 

0.92 

0.33 

0.72 

kd 

0.76 

0.22 

0.21 

0.42 

_ 

_ 

kde 

_ 

_ 

_ 

_ 

_ 

_ 

CM 

CD 

_ 

_ 

_ 

0.04 

_ 

-0.09 

Ka3 

_ 

_ 

_ 

0.03 

_ 

-0.08 

^max 

_ 

_ 

0.14 

0.06 

0.96 

1.01 

Kmb 

_ 

_ 

_ 

_ 

-0.07 

-0.08 

Kunb 

_ 

_ 

-0.10 

-0.05 

-0.81 

-0.90 

Note:  The  coefficients  less  than  |0.01|  are  not  reported. 


HED  for  Cancer  Risk  Estimates 

Under  steady  state  conditions,  the  estimated  daily  plasma  and  liver  AUC  values  in  mice 
treated  with  2. 1  mg/kg/day  of  DCA  in  drinking  water  were  both  0.05  mg/L/day.  The  average 
daily  amount  of  DCA  metabolized  per  kg  of  liver  (Ami)  was  38.5  mg/kg  liver.  In  the  present 
human  model,  assuming  a  70  kg  person  ingests  2  L  of  water  per  day  (1  L  every  5  minutes  at  12 
hour  increments),  oral  ingestion  of  0.1  mg/kg/day  (3.5  mg  DCA/L  drinking  water)  results  in  a 
free  plasma  AUC  and  a  liver  AUC  equal  to  0.05  mg/L/day.  Oral  ingestion  of  0.02  mg/kg/day 
(0.7  mg  DCA/L  drinking  water)  results  in  a  total  plasma  AUC  equal  to  0.05  mg/L/day,  while  oral 
ingestion  of  1.0  mg/kg/day  (35  mg  DCA/  L  drinking  water)  results  in  an  Ami  equal  to  38.5  mg/kg 
liver/day.  HED  values  and  the  percent  remaining  GSTzeta  activity  (PVmax)  are  shown  in  Table 
6. 
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Table  6.  Calculated  human  equivalent  doses  (HEDs)  for  different  dosimetries  of  DC  A 
_ carcinogenicity _ 


Dosemetrics  of  DCA  Carcinogenicity 

Total  Plasma 
AUC 
(mg/L/d) 

Free  plasma 
AUC 
(mg/L/d) 

Liver  AUC 
(mg/L/d) 

Amount  of 
metabolized  DCA 
(mg/kg  liver/d) 

HEDs  (mg/kg/day) 

0.02 

0.1 

0.1 

1.0 

%  Remaining 
GSTzeta  (PV^ax) 

94.4 

78.0 

78.0 

27.5 

DISCUSSION 

In  the  present  study,  PBPK  modeling  was  used  to  quantitatively  describe  IV  and  oral  intake  of 
DCA  in  humans.  The  range  of  administered  doses  of  DCA  reported  in  the  pharmacokinetic 
literature  was  exceptionally  wide  (pg/kg  to  mg/kg)  and  included  both  single  and  repeated  dose 
administration.  Collectively,  these  published  clinical  pharmacokinetic  studies  provide  an 
excellent  database  for  the  development  of  a  human  DCA  model,  especially  when  coupled  with 
information  derived  from  existing  PBPK  models  for  DCA  in  laboratory  animals.  The  historical 
focus  of  research  on  DCA  metabolism  has  been  on  the  identification  of  enzymes  responsible  for 
bio-transformation  of  DCA  and  the  mechanisms  involved  in  inhibition  of  its  own  metabolism. 
Our  present  human  model  for  DCA  suggests  that,  at  low  doses  such  as  those  encountered  in  the 
environment  (pg/L  in  drinking  water),  DCA  inhibition  of  its  own  metabolism  is  probably  a 
minor  contributor  to  the  pharmacokinetic  behavior  of  DCA,  but  plasma  protein  binding  is 
suggested  to  play  an  important  role.  For  higher  doses  of  DCA,  such  as  therapeutic  doses, 
inhibition  of  its  own  metabolism  can  be  substantial,  resulting  in  slowed  clearance  of  DCA  from 
the  body.  Metabolic  recovery  from  the  DCA  insult  appears  to  be  very  slow  and  warrants  further 
investigation. 

In  our  view,  varying  degrees  of  success  were  reached  in  describing  these  diverse  and  sometimes 
variable  kinetic  data  sets  for  DCA  with  the  present  human  PBPK  model.  The  model  parameters 
were  constrained  to  single  values  to  represent  the  mean  kinetic  behavior  for  groups  of 
individuals,  thus  the  goal  was  not  to  fit  each  data  set,  but  to  determine  model  parameters  that 
could  adequately  describe  most  of  the  data  sets.  When  inter-individual  variability  was  evaluated 
for  metabolism  and  oral  uptake  of  DCA,  a  fairly  wide  range  of  metabolic  capacity  values  were 
estimated  (Table  3).  This  variability  in  estimated  metabolic  capacity  across  individuals  may 
reflect  polymorphisms  in  GSTzeta,  the  primary  hepatic  enzyme  responsible  for  metabolism  of 
DCA.  However,  other  factors  such  as  oral  uptake  kinetics  may  be  important  as  well. 

In  the  present  simulations,  the  breakdown  of  DCA  is  assumed  to  occur  only  by  hepatic  GSTzeta 
and,  depending  on  dose,  its  metabolism  was  partially  inhibited  by  covalent  binding  of  a  reactive 
intermediate  metabolite.  This  process  is  described  as  a  second  order  process  which  depends  on 
the  concentration  of  reactive  intermediates  available  for  binding  to  GSTzeta  and  the  availability 
of  GSTzeta  in  the  liver  (Equation  4).  Experimental  evidence  demonstrates  the  dose  dependent 
inactivation  of  GSTzeta  by  DCA  (Tzeng  et  al,  2000)  and  the  subsequent  degradation  of  the 
GSTzeta  protein  (Anderson  et  al,  1999).  With  the  current  formulation  of  the  human  PBPK 
model  for  DCA,  repeated  IV  administration  of  the  lowest  reported  clinical  dose  of  DCA  (2.5 
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|j,g/kg/day)  would  result  in  less  than  1%  reduetion  in  GSTzeta  aetivity  under  steady  state 
eonditions.  With  eombined  doses  of  2  mg/kg  oral  and  0.3  mg/kg  fV  (Sehultz  and  Shangraw, 
2006),  the  model  predieted  deerease  in  metabolie  eapacity  is  nearly  13%,  whieh  persists  until  the 
seeond  eombination  dose  15  days  later,  due  to  the  slow  reeovery  of  GSTzeta  as  well  as  the  0.02 
mg/kg  DC  A  administered  daily  in  drinking  water  between  the  2  ehallenges.  If  a  person  is 
administered  50  mg/kg/day  by  IV  infusion  until  steady  state  is  reached,  only  about  13%  of  the 
initial  metabolic  capacity  is  predicted  to  remain.  Stacpoole  et  al.  (1998a)  and  Larson  and  Bull 
(1992b)  propose  that  DC  A  can  also  be  dechlorinated  to  monochloroacetic  acid;  however,  this 
alternative  pathway  has  not  been  thoroughly  investigated. 

The  mechanism  for  DCA  induced  hepatic  toxicity  and  cancer  is  unknown.  Possible  mechanisms 
associated  with  DCA  toxicity  include  competitive  inhibition  of  tyrosine  catabolism  and  the 
reaction  of  glyoxylate  (the  reactive  metabolite  of  DCA)  with  cellular  macromolecules  (Anderson 
et  al,  2004).  The  ability  of  DCA  to  cause  hepatic  cancer  in  humans  remains  equivocal. 
Interestingly,  when  DCA  was  used  as  a  therapeutic,  increases  in  hepatic  cancer  incidence  in 
humans  was  not  reported.  Also,  careful  experimental  studies  are  necessary  to  verify  the  DCA 
plasma  binding  constants  for  human  plasma.  Presently,  the  basis  for  our  model  predictions  of 
low  plasma  DCA  concentrations  is  derived  from  previous  studies  showing  low-level  binding  in 
both  rat  (Schultz  et  al,  1999)  and  human  plasma  (Chu,  1987).  However,  in  these  studies,  there 
are  not  enough  experimental  data  to  estimate  the  number  of  binding  sites  or  the  binding  affinity. 
The  binding  parameters  were  estimated  as  described  in  the  methods. 

From  a  risk  assessment  perspective,  using  model-predicted  DCA  daily  liver  AUC  as  the 
dosimetric,  the  PBPK  model-predicted  human  equivalent  dose  resulting  in  the  same  dosimetric 
as  oral  administration  of  2.1  mg/kg/day  DCA  in  mice  was  0.1  mg/kg/day.  Theoretically  this 
exposure  corresponds  to  a  10%  increase  in  incidence  in  liver  cancer  based  on  2  year  rodent 
bioassay  data  (DeAngelo  et  al,  1999).  Thus,  if  a  70  kg  person  ingested  2  L  of  water  each  day, 
this  would  equate  to  ingestion  of  3.5  mg  DCA  per  L  of  water.  The  current  maximum 
contaminant  level  (MCL)  is  60  pg/L  for  DCA  and  4  other  haloacetic  acids  (U.S.  EPA,  2002). 

In  conclusion,  the  present  DCA  PBPK  model  in  humans  is  able  to  quantitatively  describe 
the  DCA  kinetics  across  different  populations  and  for  two  exposure  scenarios.  Model 
simulations  suggest  that  environmental  exposure  to  DCA  in  drinking  water  (-0.5  pg/kg/day)  is 
expected  to  have  very  limited  effects  on  hepatic  GSTzeta  activity.  Future  studies  such  as 
genotyping  of  individual  subjects  will  allow  determination  of  the  effect  that  GSTZ 
polymorphism  has  on  in  vivo  DCA  pharmacokinetics. 
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APPENDIX  A:  DIFFERENTIAL  EQUATIONS  USED  IN  THE  PBPK  MODEL 


Nomenelature 

QC 

Cardiae  output  (L/hour) 

Cv 

DCA  eoneentration  in  the  venous  blood  (mg/L) 

Ca 

Free  DCA  eoneentration  in  the  arterial  blood  (mg/L) 

RA 

Rate  of  change  in  the  amount  of  DCA  (mg/hour) 

V 

Tissue  volume  (L) 

DIV 

iv  dose  (mg/kg) 

DO 

Oral  dose  (mg/kg) 

AO 

Amount  of  oral  absorption  (mg) 

VmaxB 

Maxim  binding  capacity  rate  (mg/hour) 

Kmb 

Affinity  constant  of  plasma  binding  (mg/L) 

CLunb 

Unbound  clearance  (L/hour) 

Kunb 

First-order  dissociate  rate  constant  (/hour) 

Kass 

Association  rate  constant  (L/mg/hour) 

Bmax 

Maximum  amount  of  binding  protein  in  the  plasma  (mg) 

RAM 

Rate  of  metabolism  in  the  liver  (mg/hour) 

cu 

Urinary  clearance  (L/hour) 

Subseripts 

pla 

_bnd 

free 

total 

1 

k 

s 

r 

GIi,  GI2 


Plasma  eompartment 

Bound  DCA  in  plasma  eompartment 

Free  DCA  in  plasma  eompartment 

Total  DCA  eoneentration  in  plasma  eompartment 

Liver 

Kidney 

Slowly  perfused  eompartment 

Rapidly  perfused  eompartment 

L*  and  GI  eompartments,  respeetively 
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The  rate  of  ehange  in  the  amount  of  free  DCA  in  the  plasma  compartment: 
dApla 

free  = 


dt 


=  QC  X  (Cv  -  Ca)  -  RApla 


bnd 


Cv  = 


Y,Q.xCv,+DIV 

QC 


i  =  slowly  perfused,  rapidly  perfused,  liver  or  kidney. 


The  rate  of  change  in  the  amount  of  DCA  bound  in  the  plasma  compartment: 

=  K  x( 

ass  V 


dt 

B^^xK. 


K^,+Ca 


CL^rb  ^  Ca 


B  X  K 

max  —mi  )  x  Cfl  “  X  Vplu  X  Cu  =  ^^22 - 22^  X  Cu  -  X  Vplu  X  Cu 


Kb  +  Ca 

where  =K^„,IK^es- 


Kb  +  Ca 


bnd 


The  total  DCA  concentration  in  plasma  compartment: 
_  Apia  f^^^+ Apia 


Ca 


total 


Vpla 


The  rate  of  change  in  the  amount  of  DCA  in  the  two  GI  compartments: 
dA 

RAan  =^  =  DO-K^,x  A^,,  -  x  A^,, 
dt 


RAgj2 


dA 


GI2 


dt 


^al  ^  Agix  X  ^G/2 


The  rate  of  oral  absorption: 

PAO  =  ^^  =  K^,xAe.„  +K^,xAe 


dt 


^a\  ''  ^^G/1  '  ''  ^^GI2 


The  rate  of  change  in  the  amount  of  DCA  in  the  liver  compartment: 
dA 

RA,=^  =  QLx{Ca-Cv,)-RAM  +  RAO 
dt 
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The  rate  of  ehange  in  the  amount  of  DC  A  in  the  kidney  compartment: 

RA^=^  =  QKx(Ca-Cv^)-CL^xCa 
dt 

The  rate  of  change  in  the  amount  of  DC  A  in  the  slowly  perfused  and  rapidly  perfused 
compartments: 
dA 

RA^=^  =  QSx{Ca-Cv^) 
dt 

dA 

RA^=^  =  QRx{Ca-Cv^) 
dt 
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APPENDIX  B:  HUMAN  DCA  PBPK  MODEL  CODE 

PROGRAM 

!  Human  DCA  Model  -  December,  2007 

INITIAL 
!  Blood  Flow  Rate 

CONSTANT  QCC  =  16.5  ICardiac  output  (1/hr-kg): Brown,  1997; 

CONSTANT  QLC  =  0.227  IFractional  blood  flow  to  liver:  Brown,  1997; 

CONSTANT  QKC  =  0.175  IFractional  blood  flow  to  kidneyiBrown,  1997; 

!  Tissue  Volume 

CONSTANT  BW  =  70  !Body  weight  (kg) 

CONSTANT  VLC  =  0.026  IFraction  liver  tissue(male):Brown,  1997; 

CONSTANT  VKC  =  0.0044  IFraction  Kidney  tissue(male):  Brown,  1997 
CONSTANT  VPC  =  0.044  IFraction  Plasma:  Merrill,  2005 

IPartition  Coefficient  for  DCA  :  mouse  values  (Abbas  &  Fisher,  1997) 

CONSTANT  PL  =  1.08  !Liver/blood  partition  coefficient:  Abbas  &  Fisher,  1997 
CONSTANT  PS  =  0.37  ! Slowly  perfused  tissue/blood  partition:  Abbas  &  Fisher,  1997 
CONSTANT  PR=  1.08  ! Richly  perfused  tissue/blood  partition:  Abbas  &  Fisher,  1997 
CONSTANT  PK  =  0.74  Ikidney  tissue/blood  partition:  Abbas  &  Fisher,  1997 

!  Metabolism  Parameters 

CONSTANT  VMAXC  =27. 16  [Maximum  velocity  of  metabolism  (mg/hr-lkg) 
CONSTANT  KM  =  6.0  IMichaelis-Menten  constant  (mg/1):  Tong,  1998 
CONSTANT  KD  =  0.3  Irate  constant  for  loss  of  enzyme 
CONSTANT  KDE  =  0.001  Ibasal  enzyme  destruction  rate(/hr):Curry,  1991 
[Proposed  in  rodent  model  (Keys,  2004).  Not  active  in  Human  Model  DCA  Model: 
[CONSTANT  KFC  =  0.0  [First  order  metabolism  rate  constant  (/hr /kg) 

[Urinary  Clearance 

CONSTANT  Clrc  =  0.0007  [1st  order  urinary  elimin.  rate  constant(/hr):Curry,  1991 
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[Molecular  Weight 

CONSTANT  MWDCA  =  129.  [Molecular  weight  DCA  (g/mol) 


['Scaled  parameters' 

[Physiological 
QC  =  QCC*BW**0.75 
QL  =  QLC*QC 
QK  =  QKC*QC 

QR  =  0.76*QC-QL-QK  [Brown,  1997 
QS  =  0.24*QC  [Brown,  1997 

VL  =  VLC*BW  [Liver  volume 

VK  =  VKC*BW  [Kidney  volume 

VS  =  0.747*BW  [Brown,  1997 
VR  =  0.129*BW-VL-VK  [Brown,  1997 
VPL  =  VPC*BW  [Plasma  Volume 

[Metabolism  and  Elimination' 

VMAX  =  VMAXC*BW**0.75 

[KF  =  KFC/BW**0.25  :  Proposed  1st  order  pathway  (Keys,  2004;  Not  active  in  Human 
Model) 

Clr  =  Clrc*BW 

[resynthesis  rate  (mg/hr)/hr  =kde*vmax  so  that  vmaxt=vmax  at  steady-state 
KS  =  KDE*VMAX 


['Dosing' 

CONSTANT 

CONSTANT 

(5  or  15  d). 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 


PDOSE  =  0.0  [ first  oral  dose(mg/kg) 

PDOSE14=  0.02  [dose  for  daily  dose  (mg/kg):  Schultz  &  Shangraw  (14  d);  Jia 


TOR  =  0.0  [  Start  time  of  first  oral  dose  (hr) 

TOR14  =  24  [start  time  of  14  days,  0.02  mg/kg  DCA 
IVDOSEl  =  0.0  [First  IV  dose(mg/kg):  For  Fox  (1996)  Data 
IVDOSE2  =  0.0  [Second  IV  dose  (mg/kg):  For  Fox  (1996)  Data 
TINF  =  0.5  [Fength  of  IV  mfusion(hrs) 

TIVI  =  0.0  [Start  time  of  fist  iv  dosing  (hr) 

TIV2  =  8.0  [Start  time  of  second  iv  dosing  (hr) 

IVINTl  =360  [Interval  of  1st  iv  repeat  dosing  (hrs) 

IVINT2  =  10000  [Interval  of  2nd  iv  repeat  dosing  (hrs) 

FinallV  =  9  [Fast  IV  dosing  time  (hr) 

FinalOR  =  337  [last  oral  dosing  time  (hr) 

ORINT  =  360  [Interval  of  oral  dosing  (hrs) 

ORINT14  =24  [Interval  of  daily  oral  dosing  (hrs) 

TGAV  =  0.05  [Fength  of  oral  gavage  mfusion(hrs) 

K1  =0.001  [Oral  absorption  transfer  rate  GIl  to  liver  (/hr) 

K2  =  3  [Oral  absorption  transfer  rate  GIl  to  GI2  (/hr) 
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CONSTANT  K3=5 


!Oral  absorption  transfer  rate  GI2  to  liver  (/hr) 


!  Plasma  Protein  binding  constants  (Clewell,  2003) 

CONSTANT  Vbp  =  91  [Maximum  Binding  Capacity  Rate  (mg/L) 

CONSTANT  CLunbp  =  9  [Unbound  clearance  (L/hr) 

CONSTANT  Kmbp  =  22.6  [Binding  Affinity  (mg/L) 

[Derived  Equations  for  Reported  Bmax  and  Kunb  (Communication  with  Dr.  Andersen  on  May 
11,2007): 

Kunb  =  CLunbp/VPL  [Dissociated  Constant  (/hr) 

Bmax  =  Vbp/Kunb  [Maximum  Binding  Capacity  (mg) 

END  [  INITIAL 

DYNAMIC 

ALGORITHM  lALG  =  2 
NSTEPS  NSTP  =  I 
MAXTERVAL  MAXT  =  L0e9 
MINTERVAL  MINT  =  LOe-9 
CINTERVAL  CINT  =  0.0 1 
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DERIVATIVE 


! - IV  Multiple  Doses - ! 

!IV  =  Intravenous  infusion  rate(mg/hr) 

!IV1:  Repeated  Dosing;  IV1+  IV2:  Fox  (1996)-  30+15,  60+30  and  100+50 
iflagl  =  PULSE(TIVl,ivintl,tinf)*PULSE(0,tstop,linaliv) 
iflag2  =  PElESE(TIV2,ivint2,tinf) 

IVl  =  IVDOSEl*BW/TINE*iflagl  I'rate  iv  dosing(mg/hr) 

AIVl  =  INTEG(IV1,0.0)  !  Amount  iv  dosed(mg) 

IV2  =  IVDOSE2*BW/TINF*iflag2 
AIV2  =  INTEG(IV2,0.0) 

! - oral  bolus - ! 

lODOSE  =  oral  gavage  infusion  rate(mg/hr)  for  lst&  2nd  ehallenge 

oflag  =  PElESE(TOR,orint,tgav) 

ODOSE  =  PDOSE*BW/tgav*oflag  I'rate  oral  gavage  dosing(mg/hr) 

AO  =  INTEG(ODOSE,0.0) 

rst  =  -(kl*ast)-(k2*ast) 
ast  =  ao  +  integ(rst,0.0) 
rug  =  (k2*ast)-(k3*aug) 
aug  =  integ  (rug,  0.0) 
rao  =  kl*ast  +  k3*aug 
aao  =  integ(rao,  0.0) 

1140DOSE  =  oral  gavage  infusion  rate(mg/hr)  for  14  d  (Sehultz) 

!  &  Daily  dosing  in  Jia  et  al.  (2006)for  15  d 

oflagl4  =  PUESE(TOR14,orintl4,tgav)*PULSE(0,tstop,finalor) 

ODOSE14  =  PDOSE14*BW  /tgav*oflagl4  I'rate  oral  gavage  dosing(mg/hr) 

A014  =  INTEG(ODOSE14,0.0) 

rstl4  =  -(kl*astl4)-(k2*astl4) 
astl4  =  aol4  +  integ(rstl4,0.0) 
rugl4  =  (k2*astl4)-(k3*augl4) 
augl4  =  integ  (rugl4,  0.0) 
raol4  =  kl*astl4  +  k3*augl4 
aaol4  =  integ(raol4,  0.0) 

! - VMAXT  in  liver - ! 

! Isotopes  (I2C-DCA&I3C-DCA)  Competitive  Metabolism  in  liver 
CVE  =  (AL12+AL13)/(VE*PE) 

RLME  =  -kd*(Vmaxt*CVE/(Km+CVE))*(Vmaxt*Km/(Km+CVE))  !  Aetive  metabolit 
&  Eree  Enzyme  Lilly,  1998 

RCME  =  RLME+KS  -  (KDE*VMAXT) 
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VMAXT  =  INTEG(RCME,VMAX) 


Iproportion  of  original  Vmax  inhibited 

VMAXP  =  VMAXT*  lOOA/MAX 


Model  for  IV  Dosed  DC  A; - ! 

—  1)  13C-DCA  in  Schultz  &  Shangraw  (2006)  and  Jia  (2006) — ! 

— 2)  12C-DCA  in  therapeutic  dose - ! 

IBlood  13C-DCA(No  Binding) 

!  Cvl3  =  (QE*CVL13  +  QS*CVS13  +  QR*CVR13  +  QK*CVK13  +  IVl  +  IV2)/QC 
!  cal3=cvl3 

!  Plasma  Compartment  (with  Binding) 

!  CA13  =  Eree  DCA  Concentration  after  IV  (mg/E) 

CVI3  =  (QE*CVLI3  +  QS*CVSI3  +  QR*CVRI3  +  QK*CVKI3+  IVI  +IV2)/QC 
RPEI3  =  QC*(CVI3-CAI3)-  RbindI3 
APEI3  =  integ(RPEI3,  0.0) 

CAI3  =APEI3A/PE 

!CabindI3  =  Bounded  DCA  Concentration 

RbindI3  =  (Vbp*CaI3)/(CaI3+Kmbp)-CEunbp*CabindI3  !(Clewell,  2003) 

AbindI3  =  integ(RbindI3,  0.0) 

CabindI3  =  AbindI3A^PE 


!  Total  Plasma  Concentration 
Apltotl3  =  APE13+Abindl3 
Catl3  =  (APE13+Abindl3)A/PE 


!  AS  =  Amount  in  slowly  perfused  tissues  (mg)' 
RAS13  =  QS*(CA13-CVS13) 

AS13  =  INTEG(RAS  13,0.0) 

CVS13  =AS13/(VS*PS) 

CS13  =AS13A/S 

!AR  =  Amount  in  rapidly  perfused  tissues  (mg)' 
RAR13  =  QR*(CA13-CVR13) 

AR13  =  INTEG(RAR13,0.0) 

CVR13  =  AR13/(VR*PR) 

CR13  =AR13A/R 


!AL  =  Amount  in  liver  tissue  (mg)' 
RAE13  =  QE*(CA13-CVE13)  -  RAMi 
AE13  =  INTEG(RAL13,0.0) 

CVE13  =  AE13/(VE*PE) 
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CL13  =AL13A^L 
AUCL13  =  INTEG(CL13,0.0) 


!  Metabolism  in  liver 

!AMli :  GSTZ-mediated  metabolism  (mg) 
RAMli=(VMAXT*CVL13)/(KM+CVL13+CVL12) 
AMli  =  INTEG(RAM  11,0.0) 


!AM2i;  non-speeific  degradation  (mg) 
!  RAM2i=(KE*CVL13*VL) 

!  AM2i  =  INTEG(RAM2i,0) 


!AMi:  total  metabolism  amount(mg) 
RAMI  =  RAMli!  +  RAM2i 
AMI  =  AMli  !+ AM2i 


!  AK  =  Amount  in  kidney  tissue  (mg) 

RKul3  =  CLr*CA13 

AKul3  =  integ(RKu  13,0.0)  [amount  excreted  in  urme(mg) 
RAK13  =  QK*(CA13-CVK13)  -RKU13 
AK13  =  INTEG(RAK13,0.0) 

CVK13  =  AK13/(VK*PK) 

CK13  =AK13A^K 
AUCK13  =  INTEG(CK13,0.0) 


ITMASS  =  mass  balance  (mg)' 

TMASS13  =(AL13+AK13+AS13+AR13+Apltotl3)+(AMi+AKul3)  !  Amount  Dosed  + 
Amount  out  (mg) 


- Model  for  oral  dosed  DC  A  (12C-DCA) - ! 

- 12C-DCA  in  Schultz  &  Shangraw  (2006) - ! 

[Blood  12C-DCA  (without  Binding) 

[CV12  =  (QE*CVE12  +  QS*CVS12  +  QR*CVR12  +  QK*CVK12)/QC 

[CA12=  CV12 

[AUC12  =  mteg(CA12,  0.0) 

[  Plasma  Compartment  (with  Binding) 

[  CA12  =  Eree  DCA  Concentration  after  oral  (mg/E) 

CV12  =  (QE*CVE12  +  QS*CVS12  +  QR*CVR12  +  QK*CVK12)/QC 

RPE12  =  QC*(CV12-CA12)-  Rbmdl2 
APE12  =  mteg(RPE12,  0.0) 

CA12  =  APE12A^pl 

[  Bound  DCA  after  oral  bolus 
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Rbindl2  =  (Vbp*Cal2)/(Cal2+Kmbp)-CLunbp*Cabindl2 
Abindl2  =  integ(Rbindl2,  0.0) 

Cabindl2  =  Abindl2A^pl 
AUC12free  =  integ(CA12,  0.0) 

!  Total  Plasma  Concentration  (mg/L) 

Apltotl2  =  APL12+Abindl2 
Catl2  =  (APL12+Abindl2)A/pl 
AUC12  =  integ  (catl2,  0.0) 


!AS12  =  Amount  in  slowly  perfused  tissues  (mg)' 
RAS12  =  QS*(CA12-CVS12) 

AS  12  =  INTEG(RAS  12,0.0) 

CVS12  =  AS12/(VS*PS) 

CS12  =  AS12A/S 

!AR12  =  Amount  in  rapidly  perfused  tissues  (mg)' 
RAR12  =  QR*(CA12-CVR12) 

AR12  =  INTEG(RAR12,0.0) 

CVR12  =  AR12/(VR*PR) 

CR12  =  AR12A^R 


!AE12  =  Amount  in  liver  tissue  (mg)' 

RAL12  =  QE*(CA12-CVE12)  +  RAO  +  RAO  14  -  RAMw 
AE12  =  INTEG(RAL12,0.0) 

CVE12  =  AL12/(VE*PE) 

CE12  =  AL12/VE 
AUCE12  =  INTEG(CE12,0.0) 


!AMlw  =  Metabolism  in  liver  after  oral 

!AMlw  :  GSTZ-mediated  metabolism  (mg) 
RAMlw=(VMAXT*CVE12)/(KM+CVE12+CVE13) 
AMlw  =  INTEG(RAMlw,0.0) 


!AM2w:  non-speeifie  degradation  (mg) 
!  RAM2w=(KE*CVE12*VE) 

!  AM2w  =  INTEG(RAM2w,0) 


!AMw:  Total  Metabolism  Amount  after  Oral(mg) 
RAMw  =  RAMlw  !+  RAM2w 
AMw  =  AMlw  !+AM2w 


!AK12  =  Amount  in  kidney  tissue  (mg)' 
RKul2  =  Clr*CA12 
AKul2  =  mteg(RKul2,0.0) 

RAK12  =  QK*(CA12-CVK12)  -RKU12 
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